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ABSTRACT 

The object of the thesis is to investigate the signal-to 
noise ratio of multi-myoelectric processors and to correlate 
the electromyographic potential with the muscle force for 
non-isometric contraction of muscles. The electromyographic 
potential, picked-up by the surface electrode is termed as 
surface emg. The non-isometric contraction has the usual 
significance that while exerting a certain force, the muscle 
length changes. Experimentally, the relationship between 
surface emg and m.uscular force has been found to be nonlinear. 
The correlations between SEMG and muscular force can be 
used to determine muscular force by measuring surface emg 
during locomotion. The amount of emg output can be used 
as a measure of the relative tension in the muscle and 
from this assessment coiild be made how load and position 
affect myogenic spinal joint compression and musc-ular 
fatigue. The following two muscles have been studied in 
detail! "extensor carpi ulnaris“,'labductor digiti minimi" 
at different muscle contraction rates. 

The subject matter of this work has been presented 
in five chapters! Chapter 1, gives a brief intrO'duction 
to general physiology excitable tissue and bioelectric 
phenomena. The effect of stimulus to cause excitation and 
muscular contraction are discussed in brief. Finally, the 
problem is outlined and form.ulated. 
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The experimental set-up for the measurement of force 
and emg at different muscular contraction rates have been 
described in Chapter 2 . It includes surface electrodes# 
multi -El^ channels# a force recorder set-up# multi-channel 
recorder# a FFT analyser and GPIB interface enabling data 
transfer between the FFT analyser and personal computer. 

The schone for measurement of force and emg at three 
different contraction velocities have been discussed. 

Chapter 3 deals with a multi-myoelectric channel 
(multi-MEC) model for improvement of the signal-to-noise 
ratio (sNR). It is shown that SNR depends on both corre- 
lation coefficient and number of channel used. The channel 
performance analysis is carried out for correlated multi-MEC. 
A time delay technique for SNR improvement has also been 
studied. In both cases it is found that SNR of multi-MEC 
is more than SNR of signal pick-up by single channel. 

In Chapter 4# both analog and digital techniques to 
linearize integrated surface emg (ISEM3) have been discu- 
ssed. The transfer function of black box in which ISEM3 
is input and the output which approximates the muscular 
force has been foiond out using both analog and digital 
technique. 

Finally scope of further work has been discussed in 


% 


Chapter 5 



chapter » 1 


INTRODUCTION AND FORMULATION OF THE PROBLEM 


1 • 1 ELECTRICAL IMPULSES IN HUMAN BODY AND ITS TRANSMISSION 
1.1.1 Introduction 

Electromyography is concerned with electric phenomena 
occurring in the nnascular tissues and its reaction to elect- 
rical stimuli. Electromyography requires the knowledge of 
electrophysiology of living tissues# bio-electric phenomena 
associated with them and a sound knowledge of electrical ins~ 
trumentation and signal processing. 

Historically# the presence of electricity in living 
organism was detected as early as in the times of Galvani 
and Volta around 1790. By using a simple galvanometer# 
Mattsucci showed in 1838 that the muscle exterior was elect- 
rically positive to muscle interior. He also proved that the 
potential difference of the state of rest declined sharply 
during excitation. He further proved that if a second nerve 
was brought in contact with a contracted muscle# it lead to 
the contraction of the second also. It can then easily be 
inferred that the potential generated due to contraction of 
the first muscle is strong enough to stimulate the nerve in 


contact. 
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1.1.2 Resting Membrane Potential and its Origin 

A potential difference of 60-90 millivolts exists 
between the outer surface of axon and its inner protopl- 
asm. This potential difference is called the resting 
membrane potential . The surface is positive with respect 
to the protoplasm. 

The chemical composition of fluid inside the nerve 
cell is different from the flxiid outside the seraipermeable 
cell membrane. The inner fluid contains a lower concentra- 
tion of sodium and chloride ions and a higher concentration 
of potassium ions. These differences in concentration of 
ions would catise a diffusion across the membrane. The mem- 
brane does not permit a free passage of ions. The cell 
contains a high concentration of protein anions which can- 
not diffuse through the membrane into the fluid outside, 
where the concentration of protein is low. Because electric 
charges of opposite sign attract each other, the protein 
anions inside the cell attract cations which can permeate 

through the cell membrane. The monbrane is much less per- 

4 - 4 * 4 

me able to Na than to K . The tendency for K to diffuse 

out of the cell according to a difference in concentration 
gradient is coimteracted by the attraction from the protein 
anions inside the cell. As long as the Cl is in diffusi- 
ons! equilibriiam, the outward flux of k"*" is slightly larger 
than the inward fl'ux and the inward fliix of Na"** is many 


tiroes greater than the ourward flxax. Hence, there must be 
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* 4 " 4 * 

a forced movement of Na outward and K inward to balance 
the diffusion. The net effect of this different distri- 
bution of ions on the two sides of the cell membrane is a 
higher concentration of anions in the interior of axon. 

A potential difference of about 70 mV according to Nerst 
equation is present over the membrane# the interior being 
negative with respect to the exterior. A disturbance of 

the resting potential difference would occur if the Na"*^ 

* 4 * 

and K ions were allowed to move according to the diffu- 
sion gradients. An inward flow of positive ions Na"^ 
would decrease the potential and cause a hypopol arlza tion 
and an outward flow of k'*’ ions would give rise to hyper 
polarization as shown in Fig. 1.1. 

1.1.3 Action Potential and its Origin 

All axons can pass from a state of physiological 
rest to one of excitation in response to a stimulus. The 
nerve and muscle tissues# which give rise to an electrical 
impulse# travelling along the manbrane in respect to a 
stimulus# are called excitable tissues. If an electric 
current of sufficient magnitude and duration is passed 
through the living tissue it gets excited and it shows in 
the form of rapid variation of the membrane potential . Ihis 
potential is called action potential . Above a certain 
threshold combination of magnitude and duration of the 
injected current pulse# the nature of action potential 



waveform does not change as shown in Fig, 1,2. In the 
figure it is clear that the potential rises from -70 mV 
to + 30 mV and falls back to -70 mV indicating a temporal 
change in the inner surface membrane’s potential from 
being negative to positive with respect to outer surface, 
Ihe initial polarization of the membrane vanishes in the 
ascending phase of the action potential# it is called de- 
polarisation phase and the descending phase is called 
repolarization phase. The duration of action potential 
may vary from 0.2 to .5 milli-sec, in the axon, 

A precondition for the appearance of an action - 
potential is critical depolarisation of the membrane, 
Dxiring indirect stimulation the critical depolarisation 
of the membrane is affected by the acetylcholine secreted 
in the nerve endings upon arrival of a nerve impulse. The 
action potential generated in the region of the myoneural 
junction is transmitted further along whole muscle fibre. 
During direct stimulation by an electric current the 
action potential (aP) in the muscle fibre arises at the 
cathode and spreads from there throughout the fibre. The 
AP lasts for 2 to 5 milli-sec# that is to say# five to 
ten times longer than that of the motor nerve fibres 
supplying it. 
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A potential difference develops across the cell 
membrane due to presence of concentration gradient of 
ions. Action potential is generated because of the 
change in permeability of the membrane due to the sti- 
muli applied, with no stimuli# the cell membrane is 
more per ne able to potassium than to sodixom. As a 
result there is a greater flow of positively charged 
potassiiim ions from the protoplasm to the extracellu- 
lar fluid than the reverse flow of sodium cations# hence 
the outer side of the menbrane is electrically positive 
w.r.t. inner surface. 

When stimulus is applied the membrane permeability 
for the sodiiim ions becomes ten times that of potassiiam 
ions. As a result the flow of sodium ions into the cell 
exceeds the outflow of potassium ions out of the cell# thus 
reversing the sign of membrane potential. This is called 
depolarisation phase. However# the increased permeability 
to sodium ions is short lived and a restorative process 
begins soon. The sodium permeability of the membrane falls. 
This phase is called inactivation and repolarisation takes 
place, 

1.1,4 Propoqation of Action Potential 

Action potential arising in an excited cell becomes 
a stimuli for the adjoining cells and causing sympathetic 
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excitation in adjoining cells. Since the magnitude of action 
potential is larger than the threshold of stimulation# an 
undamped wave of excitation travels all along the excitable 
tissue. 

The mechanism of conduction of excitation from one 
portion of the fibre to the other portion is same as the 
generation of action potential. In both the cases when the 
depolarisation reaches a critical value (the stimulation 
threshold)# an action potential is generated. When current 
is injected into a part of the excitable tissue# the depo- 
larisation begins at the cathode# an electric cxirrent flows 
between the excited (electronegative) and resting (electro- 
positive) parts of the membrane. The speeds of conduction 
may be as hi^ as 120 metres per second. The mechanism is 
shown in Fig. 1.3. 

1.2 ACTION POTENTIAL PICKED UP BY BIPOLAR ELECTRODE 

Let h(t) be the mathematical expression of single 
motor unit (will be explained s\ibsequen tly) action poten- 
tial as shown in Fig. 1.4. A specific form that provided 
a reasonable fit as given by Shenedy Ketal [ 1 ] 

h(t) = K t C2 - at) 

where a and K are constants. 

The tenporal pattern of this action potential# as 
represented by h(t) will look different if it is picked up 
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by a bipolar electrode. Due to the attenuation in the mag- 
nitude, let 1^2 ^ picked up by the first electrode 

with respect to the ground electrode as shown in Pig. 1.5. 
Since the other electrode is at a different distance, the 
AP picked up by it with respect to same ground electrode 
is represented by h 2 (t) . The bipolar electrode is conne- 
cted to a differential amplifier, whose output is the diff- 
erence (t) - h^Ct). The AP picked up by bipolar electrode 
will have both the positive and negative peaks as shown 
in Fig . 1.6, 

1«3 THE EFFECTS OF CURRENTS ON LIVING TISSUES 
1*3. 1 The Laws of Stimulation 

The electric current is an important agent which can 
give stimulus to muscle tissue if it is below a certain 
critical value above which the tissue damage occxirs. The 
effect of this stimulus is cjuickly reversible. An electric 
current may be applied with the help of extracellular or 
intracellular electrodes, the latter being more accurate 
because it prevents branching of current. 

For a stimulus to cause excitation it must have 
necessary strength, sufficient duration and steepness. 

Uie lowest strength of stimulation which can give rise to 
an action potential is called threshold of stimulation. A 
stimulus below it is called subthreshold and above it is 
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superthreshold. The threshold of stimulation is liable to 
change depending upon physiological condition of the tissue 
and method of stimulation. 

The minimum duration for which an electric current 
must be applied to a tissue to cause an excitation is in- 
versely proportional to its strength. In Fig. 1.7/ it is 
shown that a current below a minimum strength does not 
cause excitation irrespective of the duration of applica- 
tion. This minimum current is called rheobase. "nie 
time for which a cxirrent equal to rheobase must be applied 
to cause excitation is called utilisation time which 
inplies that further prolongation of the current does not 
effect the excitation. If the current is increased# the 
time needed for excitation reduces# however# this time 
imast be above a certain minimum time. It has been observed 
that rheobase slightly fluctuates depending upon different 
physiological states of rest. 

The threshold of stimulation depends on the duration 
of stimulus and the steepness of rise as well. If the 
excitation is induced by a rectangular pulse# the threshold 
of stimulation is minimum. The threshold of stimulation 
varies inversely with respect to steepness. If the steep- 
ness is less than a mininnim no action potential develops# 
no matter how great is the final current strength. This 



happens because if the rate of increase is slow# active 
changes take place inside the tissue for sufficient time 
so as to raise the threshold of stimulation. This pheno- 
menan of adaptation of excitable tissue to a slowly 
increasing stimulxas is known as accommodation as shown 
in Fig, 1,8, 

1.3.2 Effect of Direct Current on Living Tissue 

A direct current polarises the tissue. A peculiar 
characteristic of living tissue is that on application of 
direct current/ the excitation arises at the cathode when 
the circuit is made and at the anode when the circuit is 
broken. AlsO/ the threshold of excitation at the moment 
of opening the current circuit is considerably higher than 
its closing. 

The passage of electric current through a living 
tissue prodiices changes in the membrane charge. The 
region where the anode is placed/ the positive charge on 
the outside of the membrane increases causing hyperpolari— 
sation and vice-versa. In both the cases of increase 
and decrease of current/ the excitation builds up and 
then ceases exponentially, Hiis exponential drop is due 
to the capacitive behaviovir of surface membrane as shown 
in Figure 1.9, The outer and inner surfaces of the mem- 
brane form the plates and the layer of lipoids form the 
dielectric of very high strength. The pores in the 
membrane/ which allows ions to flow# makes it a leaky 





capacitx>r 


The rate of change of membrane potential depends 
upon the resistance R and capacitance C or the membrane 
time constant. 

A rise in membrane potential at the anode (passive 
hyperpolarisation) is not accompanied by change in ion 
permeability even when a strong stimulus is applied. This 
is the reason why no excitation occurs at the anode when 
a d.c. circuit is closed. On the other hand a fall in 
membrane potential near the cathode (passive depolari- 
sation) first calls es a brief increase in permeability to 
sodium ions and then a slow increase in permeability to 
potassium ions as shown in Pig. I.IO. 

The sodium permeability starts rising when the 
current reaches 50 to 80% of the threshold and rises 
further till action potential appears. The increase in 
sodium permeability does not peak immediately. First 
of all^ the depolarisation of the msnbrane at the 
cathode causes a relatively slight increase in sodium 
permeability. Then as the positively charged sodium 
ions begin to enter the protoplasm, depolarisation of 
the membrane increases, leading to consequent consider- 
able rise in sodium permeability and hence further 
depolarisation which again increases sodium permeability. 
This is called regenerative depolarisation. 



13 


Depolarisation ► 

i 

f 

Increase in sodium permeability 


Increase of ^ ^ 

sodium flow 
into cell 

It is believed that the pores through which sodium 
ions can diffuse into the cell are plugged in a state of 
rest by calcium ions (being larger in size) and these 
calcium ions pass out of the pores when the depolarisa- 
tion occxirs in response to a stimulus and thiis makes way 
for the sodium ions* 

The sodixim permeability remains increased only for 
1/lOth of a milli-sec and then it starts reducing. Uiis 
cannot be raised again by active depolarisation due to 
the property of inactivation. This inactivation of sodium 
permeability ultimately leads to triggering of the repo- 
larisation stage* 

1»3. 3 All or None Law 

This law states that a subthi«shold stimulation 
produces no excitation while threshold stimuli prodtice 
maximum excitation immediately and is unaffected by a 
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further increase in stimulus strength. The merely states 
that either no single muscle fibre is excited or all of 
than get excited together. 

1*4 MUSCLE CONTRACTION 

The striped skeletal muscles consists of fibres Vary- 
ing in length from a few millimeters to several centimetres 
and from 10 microns to lOO microns in diameter. Each fibre 
is a multinuclear structure. The fibre is enveloped in a 
transparent sheath/ the sarcolema# which appears under micro- 
scope to be structureless. 

These fibres do not exert a constant contractile force, 
but rather contract and relax repeatedly at rates as high as 
35 times/sec. The fibres are innervated in groxjips. Each 
group, which contain from 2 to 2000 fibres depending upon 
the muscle function is innervated by a single nerve axon. The 
group of fibres together with the axon and the nerve cell 
body is referred to as motor unit (MU) [ 4 ] and is considered 
the basic functional unit of a muscle. 

Contraction of a muscle can be induced experimentally 
by stimulation. Excitation of the muscle using external 
electric ctirrent stimulxis is known as direct stimulation. 
Stimulation of a Kotor nerve, leading to contraction of 
the muscle is called indirect stimulation. The extent of an 
isolated contraction of skeletal muscle varies with the 
strength of the stimulation. With threshold strength, the 



contraction is barely noticeable but increases with an 
increase in the strength of stimulxis (sufcanaxinium con- 
traction) . Above a certain strength of stimulus/ con- 
traction remains xanchanged (maximaam contraction) . 

In natural conditions a skeletal mxiscle motor 
unit usually receives from the nervous system a number 
of nerbe impulses in rapid succession. The rate at 
which these impulses axe received by the motor unit is 
called firing frequency. Under the influence of this 
rhythmic stimulation an intense and continuous contra- 
ction of muscle occxirs/ known as tetanic contraction. 

The tetanic contraction of a muscle considerably exceeds 
the maximum amplitude of an isolated contraction. Fig.l.li 
illustrates the fact that each excitation in the process 
of rhythmic stimulation causes an additional contraction 
which is suramated to the previous one, 

Ihe nerve axon transmits impulses at very great 
speed. Hence, the muscle fibres making up a notor unit 
become excited more or less at the same time. The ele- 
ctrical activity of a motor unit is illustrated by a 
serrated line in Pig, 1,12, which shows each peak corres- 
ponds to the aggregate action potential of many simxil- 
taneously excited fibres. 


Although the muscle fibres of each motor unit are 
excited synchronously in response to an efferent impulse. 
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3.1-11 Muscular contractions with 
muli of different frequency. 

Tequent stimuli cause isolated 
ntractions (I)i more frequent 
muli , incomplete tetanus(clonus) 

I, more frequent still, complete 
:anus (III) 



3.1-13 Simultaneous record of 
tetanic contraction (A) 
and action currents (B) 
of muscle 



Fig.1-12 Electrical activity of 
a single motor unit 

(A) and a whole muscle 

(B) in man 




Fig.1-14 Tetanus in an 
isolated muscle 
f i br e 
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the fibres of the different motor imits of a muscle as a 
rule function asynchronously [l,2,3,5] because they are 
innervated by various motor neurons which emit impulses 
at different frequencies and at different moments. The 
aggregate contraction of a muscle as a whole is of a 
combined character \xnder normal conditions and resembles 
complete tetanus in form even when each of the units 
function at a slow rhythm. Therefore, by asynchronous 
activity of the motor units, all movements of the body 
are smooth even with low-frequency motor impulses. The 
asynchronous activity of motor units prevent the determi- 
nation of electrical activity of an individual MU when 
picked up by xanipolar or bipolar electrodes of diameter 
3 ram. The signal picked up is a stochastic in nature when 
picked up from whole muscle as shown in Fig. 1.5, 

MU*s are having different sizes. All the motor 
units do not fire together in response to a stimuli. 

While the force is gradually increased on the body of 
the muscle more and more MU‘s come into action [l, 2, 4, 6, 7, 8, 
9,10, 11 , 12]* Large MU*s fires only if all small motor units 
have fired 13 , since different MU's contain differ- 

ent number of fibres the magnitude of AP varies. It 
has been found tout that peak to peak amplitude of wave- 
form of AP increases with the force. Ihe phenomena of 
gradual firing of MIT's of a muscle is called recruitment 
and when a MU is first recruited it fires at a fixed 



18 


frequency. For example in bicep of a man it is 7 impulses/ 
sec. A MU which is active increases its firing frequency by 
1.4 Hz/N [ 3 ] , If for example we take a force level of 
20 N# the firing frequency of this motor xmit will become 
(7 + 1,4 x20 ) 35 Hz, This will give rise to 35 Hz component 
in the emg signal. 


The number of additional motor units recruited during 
a given increment in force declined sharply at high level of 
voluntary force [2 ], The approximate relationship is 





-b K.P 

= a e 1 


Where AN = number of newly recruited motor units at 
the force level F for an incremental change of force ( F) 
and a# b are constants [ 4 ] , This suggests that even 
though the high threshold xanits generate more contraction# 
the contribution of recruitment to the increase in volun- 
tary force declines at high force levels [ 4] . 

When the force level increases# the rate of dis- 
charge of a motor unit# already active, is increased. 

This is called rate coding. Peak to peak amplitude of 
the waveform contributed to the surface emg by a motor 
unit increases near linearly with the threshold f orce [ 2 ] 
in case of isometric contraction (explain in Section 1,5,1) 

2A = K^F (1) 

where A = amplitude of action potential . F is thres- 
hold force. 



From equation (l) it is evident that there is a 
significant tendency of the units recruited with larger 
force to contribute a greater voltage to surface erag . 

If the rate coding is neglected, all the MU's of a par- 
ticular muscle will fire at the same frequency 
at all force level. Ihis will lead to a variation in 
the surface erag proportional to the square root of mus- 
cular force. But rate coding makes it linear [4, 1 4, i 5, 1 6 ] 

As the force level is increased, more and more 
synchronization among the motor xinits is obtained [17#183 
Due to the effect of synchronization, the action poten- 
tial of a particular motor unit is effected by the ele- 
ctric fields produced by all the surro\inding motor units. 
The resulting disturbances of electric potential in and 
near the muscles are very complex because the electro- 
magnetic field is continuously changing with respect to 
time and position. When plotted as a function of time, 
they constitutes the familiar interference pattern of 
electromyography as shown in Pig. 1.15 , 


1.5 SM3 - MUSCLE ACTIVITY COFUELATION 

When a muscle contracts, two parameters are easily 
measurable. One of the parameter is the force of con- 
traction and the other is the integrated electrical acti- 
vity. The first one can be measured by a suitable trans- 
ducer and the second paraneter can be processed through 
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a EK3 -channel and can be recorded. The correlation between 
these two paratoeters can# then# be obtained experimental ly 
and this result may be used fturther for better understanding 
of the physiological mechanism of muscular contraction. 

At present# it is generally agreed that the incre- 
ased force of muscular contraction is obtained by recruit- 
ing more motor units and also by increasing the firing rates 
of the motor units already active [l# 2# 4# 6# 7# 8# 9# 10# ll] • 

When a muscle contracts normally# the force of contraction 
also depends on the change of length of the muscle. Muscles 
are# elastic and hence it will be reduced in length both the 
ends are free. The equilibrium length and the resting length 
are defined as follows [1^ ] o 

Equilibritm Length 

Unstiraulated or resting skeletal muscle is normally 
xinder slight tension# since it shortens some what (20%) 
after its tendons are cut. The length of the mattached# 
relaxed muscle at which the resting tension is zero is 
the equilibrium length. 

Resting Length 

The muscle length at which maximal contraction 
tension is developed is the resting length, This length 
has been considered close to the maximal extension possible 
\ander normal conditions in the body. 



When the tension is developed in the nuiscle/ it 
may be isometric# nonisometric or isotonic condition. 

When the tension is developed but the muscle is not 
allowed to shorten the condition is called isometric con- 
traction. When the tension is developed and the muscle 
length also changes the condition is called non-isometric 
contraction. When the muscle is allowed to shorten and 
it lifts a constant weight# the contraction is called 
isotonic# since it maintains the same tension during the 
whole process of lifting. 

1.5.1 Isometric Contraction 

When force of contraction is developed in muscle# 
it can be measxired experimentally. This measured force 
will be equal to the force of contraction provided the 
contraction occxirs in isometric condition. Experimental 
evidences in case of isometric contraction were obtained 
in early fifties that the integrated electrical activity 
is having a near linear relationship with the muscular 
force [20] , Bigland and Lippold [2Q] measured the 

integrated electrical activity for the various parameters 
of the muscle. The latter was sxibjected to increase and 
decrease of length at different velocities# or the velo- 
cities were changed for the same tension or same length 
of the muscle. In every case near linear relationship 
was obtained though the slope of the lines were different. 
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Mathematical model which can justify the linear 
relationship in case of isometric contraction are avail- 
able [l«2l,22/ 23 ] , Theoretical studies by many authors 
have attempted to establish this linear relationship owing 
to its in^ortance in practical application such as the 
proportional control of the gripping force of a prosthetic 
hand or in tactile feedback. Most of these stxidies suggest 
that the amplitude of the EMG should increase as the square- 
root of the tension rather than linearly# if no tor \inits 
fire independently of one another. 

1»5.2 Non-Isoraetric Contraction 

If change in length of the muscle is allowed# a 
passive elastic tension is developed in the muscle which 
is to be subtracted from the total measured force to obtain 
the force of contraction developed in the muscle. This is 
explained with the help of Pig. 1.16 and Fig, 1.17 , When 
an iinstimulated skeletal muscle is stretched the passive 
elastic tension increases as an exponential function of 
length over a length upto 200% of the equilibrium length as 
shown in Fig, 1.16 {in the figure# lOO represents the equi- 
libriiam length# resting length is about 12 5) [l9 3. Stretches 

upto atleast 150% of ec[uilibrium length are perfectly rever- 
sible# the muscle snaps back to its equilibrium length when 
released. Skeletal muscle ruptxires at about three times its 
equilibrium length. Figure 1.17 shows the relation between 
tension developed during maximal voluntary effort and the 



Tension (lbs) "n Tension 
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IG-116 STATIC LENGTH TENSION DIAGRAM OF 
ISOLATED SKELETAL FIBRE AT REST 



FIG-1T7 ISOMETRIC LENGTH TENSION SUMMARY 
FOR HUMAN TRICEPS MUSCLE 



length of the triceps muscle in man. The net or active 
voluntary tension cuirve is obtained by stretching the 
muscle fibres and connective tissues to any given length 
and then subtracting the value obtained from the total 
tension exerted by the contracting ni\iscle at the same 
time. 

It is obvious from these two figures that the 
force of contraction of a muscle is more difficult to 
calculate if the muscle is subjected to change in length. 
For involve testing of the muscles the experimentation 
is HK>re complicated and the precise knowledge of anatomy 
is needed. 

1*6 FORMULATION OF PRESENT PROBLIM 

A general background is given in the previoxis part 
of this chapter with emphasis on recent developments. In 
light of this the following problems have been identified 
and formulated for the present work. 

It is difficult to provide a mathematical model 
correlating the surface emg with muscular force in case 
of non-isometric contraction because it is difficult to 
take into accoxint the passive elastic tension developed 
in the muscle and also to determine it experimentally. 

So, one of the objecti"e of the present work is to give a 
transfer function relationship between the integrated EMG 
as input and force measured by transducer as output. The 
detailed experimental arrangement is given in Chapter-2. 
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The EM3 picked up by sxirface electrode in the case of 
non-isometric contraction depends upon 

(1) Bie length of muscle 

(2) The choice of muscle 

(3) The location of electrodes 

(4) Temporal pattern of emg picked up by electrode. 

1*6.1 The Choice of Muscle 

On the entire hiiman body there are hardly a few 
joints which are controlled by single muscle. In most 
of the cases# several muscles are involved for sending 
joint. In such a case# the amount of force exerted by 
muscle is difficult to know without precise knowledge of 
anatomy of all the muscles involved. 

In the present analysis the muscle# “extensor carpi 
ulneris* and “abductor digiti minimi” were selected for 
study. The muscle "extensor carpi ulneris" was more or 
less solely responsible for wrist abduction. Except for 
a little contribution (about 2094) by"fluxor carpi ulneris"# 
other muscle contribute hardly anything for this movement. 
The muscle 'Abductor digiti minimi" which is a small muscle# 
responsible for the abduction of little finger was sele- 
cted. For abdxjction of more than 2 mm (more than 1 cm in 
present study) the length of the muscle is changing# so 
the condition is non-isometric construction of nniscle. 
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1.6.2 The Length of Muscle 

A small muscle is having larger nximber of small motor 
units and fewer number of large motor units e. g. eyelid 
muscle# muscle in forearm where as a big muscle has small 
number of small MU and large number of large MU e.g. back 
muscle# thigh muscle. So in the case of small mxiscle the 
incremental force is small for each newly recruited motor 
xinit. So# it is possible to get small incremental force in 
case of small muscle as shown in Fig. 1. l8. 

1*6.3 Location of Electrodes 

Surface electrodes used to pick up eng have a dia- 
meter of 3 mm. Hence# these electrodes pick up emg only 
from a small area aroxind the electrode. When independent 
output from adjacent muscles or from closely spaced areas 
in one muscle are required# spatial selectivity may be a 
problem. An finally a f^lndamental limitation of all surface 
electrode systans# regardless of design is that only super- 
ficial muscles can be used. The action potential picked up 
by electrodes depends upon the depth of muscle fibre. The 
maximum height corresponds to the action potential coming 
from the muscle fibre nearest to the skin. Similarly# the 
minimum one is coming from the muscle fibre lying in the 
deeper layer. The variation is given by 
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exp (-0.2 r ^ 

depth of muscle studied in cm 
mid -band frequency of EMG spectra 
potential difference observed at distance 
r# radially from the fibre 
potential difference at the fibre. 


When EMG interference pattern is picked up by the surface- 
electrode# the action jxjtential of motor xinits lying near 
the skin will contribute more high frequency components. 
Motor units lying at deeper layer will contribute less high 
frequency conqponents as shown in Fig, 1.19, 

1»6»4 Temporal Nature of Action Potential 

The action potential picked up by bipolar electrode 
has both positive and negative peaks# superir^xssed on it 
the noise signal [ 24 ] . The amount of noise present in 

normally processed EMG output is relatively large. This 
noise increases in magnitude with increasing contraction 
level. 

The most commonly used technique for SNR enhance- 
ment of myoelectric channel is by applying an averaging 
filter. It is clear that the averaging filter in the ME 
channel can only increase SNR at the cost of an increase 
in response time. However# in a myoelectric control 
system or other applications the time is an in^rtant 
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factor and sufficient averaging to achieve a high SNR may 
be prohibited. So, in order to enhance the performance 
signal processing time must be small. So^ the effect of 
multi myoelectric channel on SNR in^rovement is analyzed 
under the assumption that the signals from N channels are 
uncorrelated. This is under the assumption that the signal 
from one channel is distinct from other channel electrodes. 

At closer distances# these N channel signals become corre- 
lated due to the fact that the some common motor units may 
be picked up by different channels. A model to acco\int for 
the correlation between signals in the multi-channel case 
and detailed experimental results are discussed in Chapter-3. 
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CHAPTER - 2 


THE EXPERIMENTAL SET-UP 


2 .1 INTRODUCTION 

The science of recording electrical phenomena which 
accompanies muscular activity is known as electromyography. 
This can be recorded using surface electrodes or using micro- 
electrodes which can be injected into the nniscle. In order 
to extract the EM3 signals which are deeply immersed in 
noise and are relatively weak signals/ it is required to do 
some sort of S/N ratio enhancement. When a muscle contracts 
two parameters are easily measurable, the force of contract- 
ion and integrated electrical activity. 

2.2 THE EMG ELECTRODES 

The EM3 electrodes are meant to pick up the myoele- 
ctric signals which are in the microvolt-millivolt range. 
Since the signals are of low magnitudes, extra care needs 
to be exercised in the preparation of the electrodes. 
Different types of electrodes are used in elec trophy isology 
depending upon the requirements for accuracy, reproductibility 
of results etc. Generally metal electrodes like wires and 
plates of silver# platinum, nickel etc. are used. These 
electrodes have low resistances. But it is not recommended 
to use them for long stimulations since ions from the 
electrodes pass into the tissue under the action of current 
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and may have toxic effect. Also while recording d.c. poten- 
tial# there is a chance of polarisation of these electrodes 
due to foirmation of electrochemical cellso 

Silver was chosen as the electrode materials since it 
offers the best performance to cost ratio. Silver with an 
outer layer of silver chloride is best for electrode mater- 
ial. Each active electrode was 3 mm in diameter while the 
neutral electrode was 10 ram in diameter. The silver electrode 
buttons of the two active electrodes were fixed on a plastic 
base with 1 cm separation .The neutral electrode was fixed 
separately on a plastic base. This was done to avoid corirup- 
tion of the signal by low frequency noise due to movement of 
electrodes while doing physical task. Coaxial cables of 
smallest available core diameter were chosen to transmit 
the signals picked up by the electrodes in order to minimise 
the external noise entering into the system. 


2.3 THE EM3 CHANNEL 

The signals picked up from the surface of the skin 
have to be processed in several stages before operating a 
control unit. The maxiraxm r.m. s. voltage obtained between 
two skin electrode plate of about 1 cm separation is between 
50 to 1000 UV. To obtain an output d.c. signal within a 
dynamic range of a few volts# the required gain from r.m.s. 
to d.c, mxxst be of the order of 60 to 90 ,db. The noise 
signals appearing at both electrodes under identical phase 



and amplitxide conditions must be suppressed efficiently in 
order to reduce disturbances. The suppression is obtained 
by a high CMRR and a high input inpedance in the first 
differential stage. The CMRR should be more than 70 db 
and input impedance should be more than 100 mega-Ohms. 

The interface impedance between tissue and electrodes may 
fliictuate quite considerably. If the input impedance of 
the amplifier is not high compared to the interface impe- 
dance, the fluctuations will cause troublesome errors when 
common mode signals occxir. 

The EM3 channel is primarily used to differentially 
amplify the input signals and filter the frequency regions 
of interest. The EMG channel used consisted of 

(1) the differential amplifier 

(2) the band pass filter 

(3) the adjustable gain amplifier 

(4 ) the half wave rectifier 

(5) the full wave rectifier 

(6) the low pass filter. 

The block diagram of the EM5 channel is shown in 
Fig, 2.1 and its electronic circuit realization is shown 
in Fig. 2.2. 

The differential amplifier was incorporated in the 
circuit to give noise free signals. Ihe high CMRR was 
accomplished utilizing a potentiometer in the differential 
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stage. To avoid problems arising from polarisation voltages 
and the late a.c. coupling, the voltage gain of the differ- 
ential stage has been kept low. Ihe differential amplifier 
used three inputs from the three surface electrodes (one 
reference and two active ones). It amplified only the 
signals from two active electrodes based on the potential 
difference between each of thsn and the reference electrode. 
Since any interference from external causes produced identi- 
cal changes in the potentials of all the three electrodes# 
the output signal was xinaffected by external noise. 

In the following two stages# the myoelectric signal was 
amplified and passed to a band pass filter. A resonant 
type band pass filter was used# lower cut-off frequency 
being 110 Hz and higher cut-off 280 Hz. 

The variable gain amplifier was incorporated to vary the 
gain to desired level to prevent saturation of the signal. 

The half wave rectifier was incorporated to clip off 
the negative portion of signal. The half wave rectified 
output was used to analyse the frequency characteristics of 
the EI^ spectrum. 

The full wave rectifier and subsequent low pass filter 
havii^ time constant of 100 milli-sec. was used to reduce 
ripple to a negligible value. 
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2.4 THE FORCE RECORDER 

The arrangement for the measurement of force is shown 
in Fig. 2.3. Two strain gauges (300 Ohms# ) 

used in the Wheatstone bridge configuration, as shown in 
Fig. 2.4. were mounted on two sides of a thin steel -strip. 
Since the mass of the steel -strip was small, its frequency 
response is high which permits the measurement of force at 
high muscular contraction rates. The present study is on 
non-isoraetric contraction of muscle and hence the bending 
of the strip was limited upto 0.8 to 1 cm. For this con- 
figuration of Wheatstone bridge, the sensitivity of the 
system gets doubled. This is due to the fact that the 
bending of steel strip leads to increase in resistance of 
strain gauge mounted on one side of steel strip and decirease 
in resistance of strain gauge mo-unted on other side. This 
conf igiiration is also insensitive to temperature changes 
because strain gauges used are identical. 

For the configuration shown in Pig, 2.5, even though 
it is insensitive to temperature change, the sensitivity 
is half of the conf igxiration as shown in Pig. 2.4, 

The output of the Wheatstone bridge was made to pass 
through differential amplifier to remove the effect of noise 
signal present in the mechanical part of the setup. The 
zero adjustment of the bridge was achieved with the help 
of differential amplifier and summing amplifier as shown 
in Fig, 2,4. when the wheatstone bridge is unbalanced 



Constant current source 



FlG.2-4 CIRCUIT DIAGRAM OF FORCE RECORDER 
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, Frequency response 

Tape speed 



. frsqueinCY response of the multi- 
table 2.1 ♦ channel recorder at different 

tape speeds 



fig .-2,5 
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due to bending of steel strip an output voltage appears. 
Constant current source excitation was used for the wheat- 
stone bridge so that output unbalanced voltage is also d.c. 
in nature. 

The maximxim bending of the steel strip was limited 
upto 0.8 to 1 cm because the semiconductor strain gauges 
shows linearity upto that bending limit only. Fig. 2.6 
shows the characteristic. 

The force recorder was calibrated by applying known 
forces and output voltages were measured. The calibration 
curve is drawn# as shown in Fig. 2.7, The characteristic 
in Fig. 2.7 shows that good linearity is obtained upto 
1000 9^ force. 


2.5 THE MULTI -CHAI^TNEL RECORDER 

The multi channel recorder used in the present set-up 
was a Tech - MR - 40 model which can record and reproduce 
4 channel analog signals on cassette tape. iSie principle 
features of this recorder are: 

1. 4 tape speeds# 4.8 cm/s# 2.4 cm/s# 1.9 cir/s# 1.2 cm/s 

2. channel monitoring 

3. an accurate electronic tape counter with LED display. 

4. identification code recording# reproduction and 
search. 

The main specifications of the multichannel recorder are 




FIG- 2-6 STRAIN GAUGE CHARECTERI5TIC 
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given below* 

Input impedance - 1 M Ohm^ 
Output impedance - 75 Ohm 
Output current - 10 mA., 


The frequency response of the multichannel recorder at 
different tape speed is shown in Table 2,1. A tape speed of 
4.76 cm/sec. was chosen since the frequency of response was 
sufficient for truthfully recording the EMO signal. 

2.6 THE FFT ANALYSER 

The AD-3521 is a two channel FFT spectrum analyser 
capable of very high accuracy. It has a frequency range 
of d.c. to 100 KHz# a 15 bit A/D convertor and a full 
90 db dynamic range. Both input channels are isolated to 
allow two signals of differing voltage levels to be analy- 
sed at the same time. The following fxanctions are presents 


Spectrijro 
Histogram 
Time-waveform 
Au to -c o rr el a tion 
Cross-correlation 


Transfer function 
ImpxiLse response 
Coherence function 
Coherence output power 
Cros s-spectrum. 


The AD-3521 FFT analyser has a basic memory of 2 50 Kb, 

This can be increased to 2M bytes# by the addition of option 
boards# allowing long term transient capture for future analysis 
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single or dual 5 1/4" floppy disk drives can be connected# 
as an additional option# giving storage facilities for 64 
panel conditions and 64 displays per disk. 


The main specifications of FFT analyser is given below 


1 ) 

2 ) 

3) 

4) 

5) 


6 ) 


Sampling points 
Sanpling frequency 
Frequency range 
Windows 

Delay between 
channels 


Frame node 
single frame 

Multi frame 


2048/frame 

4,096 times frequency range 
1 Hz to 100 KHz 

Rectangular# Hanning# Minimxm 
Flat top 

Single frame node-delay within 
1 frame 

Multiframe node - delay within 
the set nxjmber of frames. Delay 
is set in units based on seconds 

One frame (2048 time data item) 
is input and processed 
Multi frame of data are input. 
One arbitrarily selected single 
frame is successively extracted 
for processing. 


The analyser can also be connected to an external 
computer and print-out facilities can be obtained by conne- 
cting a video printer to the appropriate output connector 
on the rear panel of the instrument. 



2.7 THE SCHEME OF EXPERIMENT 


The scheme is shown in Fig, 2.8, Bipolar sxirface 
electrodes were used to pick up emg signal. Since one 
of the objective of the present experiment was to study 
the signal to noise ratio improvement using multimyo- 
electric channels# so two electrode pairs were motinted 
for study of “abductor digiti minimi" responsible for 
finger abduction and three electrode pair for "extenser 
carpi xilnaris". The emg signals picked up by bipolar 
electrodes have both positive and negative peaks and 
superimposed noise. The emg signal picked up by all 
different electrode pairs were summed up using a summing 
anplifier. The summed emg and emg picked up by different 
electrode pairs were recorded for study. 

To study the surface emg and muscle force correla- 
tion in case of non-isometric contraction# the sximroed 
emg was full wave rectified and low pass filtered to get 
d.c. voltage. 

Force was exerted on the lever on which strain-gaxiges 
were mounted. The deflection was limited to one cmn. 
for non-isometric contraction of the muscle. With the 
change of frequency in function generator the first light- 
spot on the oscilloscope traced saw-tooth wave with 
different velocities. Force was exerted on the lever in 


such a way that the second light spot closely follow the 
first. In this way surface en^ and force were recorded 
for different nniscular contraction rates. 









CHAPTER - 3 


A SIGNAL~TO~NOISE INVESTIGATION OF A CORRELATED MULTI 
MYOELECTRIC CHANNEL 


3,1 INTRODUCTION 

When surface detected electromyographic activity is 
to be used for some purpose such as controlling a powered 
orthotic/prosthetic device# it is often desirable that 
the control signal derived from the EMG activity be as 
smooth and responsive to contraction level changes as 
possible. Surface-detected EM3 closely resembles ampli- 
tude-modulated random noise and it is necessary that 
some initial processing be performed to convert this raw 
EMG activity into a useful signal. Ideal output charac- 
teristics of the processed EM3 would include a signal 
monotonically related to instantaneous muscle contraction 
with no superimposed noise, 

A major source of difficulty in proportional EMG 
control at present is due to relatively large amount 
of noise present in normally processed EM3 output. One 
desires a processed signal which is a smooth image of the 
modulating envelope for purpose of correlating electrical 
and mechanical events. This smoothing can be achieved by 
low-pass filtering but only at the cost of distortion in 
the amplitude and phase of modulating envelope which leads 
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to sluggish system response detrimental to control. Control 
of sluggish systems is inherently difficult for human 
because of the temporal lag and distortion in feedback 
between the user's actuation and system response. Thus 
the problem is how to obtain a smooth measure of the overall 
electrical activity of a muscle during both slow and rapid 
movements without frequency-dependent filter distortion. 

There are several methods by which the SNR can be incre- 
ased. 

Method is One approach to solving this problem is to use 
spatial averaging of several uncorrelated EM3 signals from 
a muscle, instead of temporal averaging of a single signal. 

Method 2s Using two channels and giving a time delay to 
signal of one channel and then summing. 


3.2 correlated MULTIMyOELSCTRIC CHANNEL MODEL 
3. 2.1 Model 

When multimyo electric channel is used for SNR improve- 
ment, it is ass\Jmed that the signals from N channels are 
uncorrelated. This is an approximation of the case where 
the electrode used for obtaining the signal for one channel 
is distant from other channel electrodes. At closer dis- 
tances# these N channel signals become correlated due to the 
fact that the some common motor xinits may be picked up by 
different channels. A model to accoxint for the correlation 
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between signals in the multi-channel case is discussed 
below. 

The multi - MEC model that accounts for the corre- 
lation between channels is shown in Fig, 3,1, 

Let Ujj^(t) (j = and i = 1,2....N) repre- 

sent the jth motor \anits innervation signal in the ith 
channel. Yj^Ct) represent the ith myoelectric signal. 

The mean ECy^^Ct)) of each channel output is denoted by 
and var !yj^(t)} by 

N 

Z = 2 Yi 

i=l 

N jg 

E [ zCt) ] = E [2 Yi ] = 2 i (1) 

i=l i=l 

E [2^(t) ] = e[ ( 2 y.Ctr ] 

i=l 

N 2 N 

= 2 E( y. (t)) + 22 E(yi(t) YvCt)) 

i=l i/K ^ 

( 2 ) 

E !yi(t) ] = (3) 


Cross covariance 


^ik = 

Substituting eqns. 
E [ 2^ (t) ] * 


E [ y^it) y)c(t) ] - (4) 


(3) and (4) in eq. (2) 

N 2 2 N 

2 ( + 1^4 ) + 22 ^^-5 V + M'j ^Tf) 

i=l ^ i^ ^ ^ 

C5) 
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where is covariance between channels i and k. 

Cross correlation coefficient 


^ik 

If 

Q 

O 

Q H- 


(6) 

From eg. (5) 

and eg. (6), 



E [ (t)] = 

N 2 

S OT + 

i=l 

N 

2 

i=l 

2 N 

+ 22 o 

±/^k ^ 



H 



+ 

2 

i/k 



N 2 

2 0 + 
i=l 

N 

22 

i/^k 





(7) 

1 

R^Co) - 




E [x^(t) ] - 

>^1 


R^(o) is autocorrelation at zero. 

Prom eg. (7), 




Var [ Z(t)] 

N 2 


N 

°^i \ ^ik 

Mk 1 IK 




(8) 

Frcsn eg. (l) 

and (8) 
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) 
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3.1 * A Corrolatod MulU-myootoctrIc Channel Model, where 
ond lUuli) are the common units In dllforcnt channels. 
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where SNR^ denotes the signal to noise ratio of the corre- 
lated multi MEC, 

3.2.2 The Experiments 

The experiment was conducted on two muscle one small 
"abductor digit! minimi" and other medium sized muscle 
''extensor carpi ulnaris. Two channels were used in case of 
former because it was difficxilt to place more than two 
electrodes pairs with sufficient spacing so that emg pick- 
ed up by the electrode pairs have minimum correlation. 

Ihree electrode pairs were placed in case of "extensor 
carpi ulnar is''with a spacing of 1 cm between two electrode 
pairs and the ground electrode was placed 3 cm away from the 
active electrode pairs. 

The scheme for recording surface emg picked up by 
individual electrode pair and their sum for various contra- 
ction rates of muscle has been explained in Chapter - 2 under 
“The scheme of experiment." 

The surface emgs recorded were rectified to get 
substantial mean value of the SEMG and were stored in the 
memory frame of FFT analyser. The frame to frame analyser 
was carried out using FFT analyser and various parameters 
mean, autocorrelation/ cross correlation between different 
channels were calculated. Tie rectified SEM3 picked up by 
each individual electrode pair, their spectruros# autocorre- 
lations/ cross correlation for "extensor carpi ulnaris" are 
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FIG. 3.3 : RECTIFIED EMG PICKED UP BY ELECTRODE PAIR B 
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MUSCLE CONTRACTION RATE OF 2 N/sec. 
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shown in Figs- 3-2 to Fig. 3.8, 

3.2.3 R esults and Discussion 

The results obtained from the analysis of the experi— 
data 

mental^for SEMG picked up from "abductor digit! minimi" at 
different contraction rates are shown in Tables 3.1 to 3. 3. The 
variation of SNR of surface emg picked up by each indivi- 
dual electrode and the summed signal with memory frame 
number are shown in Fig, 3.9 to Fig. 3.18. It is evident 
from the plot that the SNR of sxjmmed signal is higher than 
the SNR of SEM3 picked up by each electrode pair. Also 
the SNR increases with the memory frame n\miber which is due 
to the fact that the amount of signal present increases 
with the increasing contraction {force applied). Ihe 
Variation of SNR with muscular contraction rate is shown 
in Pig, 3.19 and Fig. 3.20. As the muscle contraction 
rate is increasing the SNR also increases because with 
increasing rate of doing work the amount of sEMG present 
in signal picked up by electrodes increases. 

3.3 SNR IMPROVEMENT BY DELAYED SUMMATION 
3.3.1 I n tro duct ion 

The most important conponent of noise present in 
surface EH5 has frequency of 5o Hz (power line frequency). 

An improvement in signal to noise ratio can be achieved 
by time delaying the surface emg picked up by one electrode 
pair and then summing it with the surface emg picked up by 


TABLE 3.1 * Signal to Noise Ratio of SEMG Picked-up from "ABDUCTOR DIGITI MINIMI" for 

Contraction Rate 0.5 N/sec. 
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TABLE 3.3 t Signal to Noise Ratio of SEM3 Picked-up from “ABDUCTOR DIGITI MINIMI 
for contraction Rate 2N/sec. 
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EJG. 3.14 1 SIGNAL TO NOISE RATIO OF "EXTENSOR CARPI XJLKARIS 

FOR CONTRACTION RATE 1 N/sec. 
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other pair of electrode. The amoiint of delay to be given 
is 10 msec# this is explained as follows. Ihe time period 
of 50 Hz noise component present in SEM3 is 20 msec. When 
a delay of 10 msec, is given to the SEM3 picked up by one 
electrode pair and is added to SEM3 picked up by other ele- 
ctrode pair# the positive peak of noise of one channel false on 
negative peak of other channel cancelling the power frequency 
noise. 

3.3.2 The Experiment 

The delay was introduced to SEMG pick up of one 
channel by FFT analyser and both the channel EMGs were 
stored in memory of FFT analyser. Ihe analysis was 
carried out on SEMGs picked up at two different contraction 
rates as shown in Pigs. 3.2l to 3.28. The data stored in 
the memory of FFT analyser were transferred through GPIB 
interface to IBM. PC, 

3.3.3 Resultsand Discussion 

The cross correlation coefficient and SNR of the 
delayed and summed SEMG has been compared with that of 
without delay and summed SEM3 and the results are shown in 
Tables 3.4 and 3, 5, It is evident that there is considerable 

improvement in the SNR in the case of delayed and s^lmmed 
compared to without delay and summed signal. There is also 
considerable decrease in cross-correlation coefficient. When 
the Fourier transform of the delayed and siimmed signal was 
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707 B'f ELECTPOOE PRlP^m CONTPRCTiON PftTE ] N/sec MISC MENU 

E-3 ’ FUNCTION 



FIG. 3.21 i EMG PICKED UP BY ELECTRODE PAIR A FROM 
MUSCLE "ABDUCTOR DIGIT! MINIMI" FOR CON- 
TRACTION rate 1 N/sec. 
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t EM3 PICKED UP BY ELECTRODE PAIR B FROM 
MUSCLE "ABDUCTOR DIGITI MINIMI" FOR CON- 
TRACTION RATE 1 N/sec AND DELAY OF 10 milli- 
seco 


FIG. 3.22 
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FIG. 3.24 : SUM OF EM3'S PICKED UP BY ELECTRODE PAIRS 

A AND B WITHOUT DELAY TO CHANNEL B 
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FIG. 3.25 : EMG PICKED UP BY ELECTRODE PAIR A FROM 

MUSCLE "EXTENSOR CARPI ULNARIS FOR CON- 
TRACTION RATE 4 N/sec. 
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FIG. 3.2 7 ; SUM OF EMG'S PICKED UP BY ELECTRODE PAIRS A 
AMD B WITH CHANNEL B DELAY OF 10 milli- 
sec. 
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taken, it was foxmd that the magnitude of 50 Hz frequency 
component also reduces. The results are shown in Table 
3.4 and Table 3.5. 

The advantage of this method in SNR improvement is 
that the SNR which can be obtained using multin^o elec trie 
channel for a given cross correlation coefficient, the 
same value of SNR can be obtained using two channels and 
delay circuit. This reduces the difficulty in placing 
multiple electrode pairs on small muscle and also reduces 
the circuit complexity. 

In the present experiment using two channels and 
delay of 10 msec, the SNR and cross correlation coefficient 
were 1.172 and 0,04 respectively for muscle contraction 
rate of 2 N/sec. For the same value of SNR and using 
multimyoel ectric channel it requires 6 channels {when 
compared with Fig. 3.29, resxol 1:3 obtained by Zhang, Parker, 
Scott [28] using multimyoel ectric channel). 



table 3.4 * Signal to Noise Ratio of SEMG Picked-up from "ABDUCTOR DIGITI MINIMI 

for Contraction Rate IN/sec. 
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TABLE 3.5 * Signal to Noise Ratio of SEMG Picked-mp from "EXTENSOR CARPI ULNARIS" FOR 

Contraction Rate 4 N/sec. 
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CHAPTER - 4 


SURFACE EMG AND FORCE CORRELATION FOR NON-ISOMETRIC 

CONTRACTION 


4.1 INTRODUCTION 

In prehistoric period cave man used to carry boulders 
in front (as seen in wall paintings) but present day field 
workers are seen carrying load on their back also. In 
order to find out the stress-strain distribution on hximan 
skeleton and the most suitable posture of carrying load, a 
suitable device to measure force during dynamic condition 
is necessary. Surface emg is one of the easiest way of 
measuring muscular force during locomotion. The peak to 
peak amplitude of action potential increases with increas- 
ing force. It has been explained in Chapter - 1 that inte- 
grated electrical activity is having a near linear relation- 
ship with muscular force in case of isometric contraction. 
Mathematical models explainir^ this near linear relation- 
ship is available owning to its importance in practical 
application such as proportional control of the gripping 
force of a prosthetic hand, dentistry etc. 

Under dynamic condition the muscle length changes 
and the contraction is non-isometric . A passive elastic 
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tension is developed in the muscle which is to be siibtra- 
cted from the total measured force to obtain the force of 
contraction developed in the muscle. It is difficult to 
measure this passive elastic tension developed in the 
muscle* So# instead of giving a mathematical model showing 
the effect of firing frequency# recruitment and rate codino 
(explained in Chapter « i) on force of voluntary contraction 
in case of non-isometric contraction# an attempt was made to 
determine muscular force from directly measuring isEMG. 

For determining force from ISEM3 a suitable black box 
has to be placed between ISEMG as input and force as output. 
Ihe transfer function of this black box is to be determined. 
When isEMG picked up during dynamic condition is passed 
through this black box# the resulting output gives the 
measure of muscular force. The block diagram of the model 
is shown in Fig, 4.lo 

The realization of the transfer fxmction can be done 
using both analog and digital technique, 

4.2 THE EXPERIMENT 

The experiment was carried on two different muscle# 
one small "abductor digit! minimi" and other medium sized 
"extensor carpi ulnaris". The reason for selection of 
these muscles has been explained in Chapter - 1. The 
force exerted by muscle and ISEM3 were recorded at diff- 
erent muscle contraction rates. 



78 


For the realization of transfer fvmction of the black 
box# the ISEMG and force recorded at different muscular 
contraction rates of two different muscle were studied. 

The ISEMG and force exerted by muscle for two different 
contraction rates are shown in Figs, 4, 15,4.19,4. 23 and 4,24, 

The autocorrelation of force and ISEMG are shown in 
Fig. 4.5 and Pig, 4,6. Ihe cross correlation between ISEMG 
and force exerted by muscle is shown in Pig., 4,7, 

In almost all cases it was found that the isEMG satur- 
rates for gradually increasing force. Fig, 4.8 and 
Fig, 4, 11 {the force is increasing linearly). This has 
been explained with suitable mathematical model for iso- 
metric contraction. The ISEMG is to be linearized 

so that it approximately represent the muscular force. 


4.3 ANALOG CIRCUIT REALIZATION 

The plot of force Vs. ISEMG is shown in Fig, 4,2, 

The mathematical expression which is found to fit the 
curve is given by 

F = exp (3.833 Ejj, ) * 0.2226 

The realization of this force Vs, ISEMG relation 
was done by antilog amplifier in which ISEMG serves as 
input and the resulting output which is linear, represents 


muscular force 



The block diagram of antilog amplifier is shown in 
Fig. 4.3 and the circuit realization is shown in Fig, 4, 4, 

and are operational amplifier with infinite 
input resistance and zero differential input voltage 


V2 = -Vf + 


Vy (In - In Iq) + 


A,.,. , — , Y 

Rl + R2 s 


( 1 ) 


V 2 is the negative of the voltage across D^. 
V 2 = - (In I 2 - In I^) 


( 2 ) 


Combining eq. (l) and eq. (2) 


R. 


S Rx + R 2 


Vj In 


Since V = I, R* 
o ^ 


from eqn, (3) it follows that 


In 


If R' 


V_ 


(3) 


= R‘ I, ln“^ - V„ ( 




1 ^ 


) (4) 


v„ 


Comparing eqn, (4) with the force-lSEMG relation# values 
of 

= 300 milli-anp# R. = 3,3K # R- =5. 


R* 


33 Ohm. 


since the input voltage to the antilog amplifier is 






Constant current source 
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limited between O to 100 roV# the ISEMG was attenuated and 
given as input. ISEMG recorded at two different muscular 
contraction rates were given as input to antilog amplifier# 
Fig. 4.8 and Fig. 4,11 and the resulting outputs are shown 
in Fig. 4.9 and Fig, 4,12. It is seen that ISEMG gets 
linearized. When the output of antilog amplifier (Fig.4.9) 
for a given ISEMG as input (Fig. 4.8) is compared with the 
force recorded (Fig, 4.10) during experiment using force 
recorder, it is seen that they are of similar nature. 

4 . 4 DIGITAL TECTNI QUE FOR THE REALIZATION OF TRANSFE R 

FUNCTION 

The transfer function of black box in which ISEMG 
serves as input and force is output is shown in Fig, 4, 13 
for muscle contraction rate of 4 N/sec . The transfer 
function was found using FPT analyser. Ihe approximate 
nature of transfer function is shown with dotted line in 
Fig. 4.13. The general nature of transfer function is 
shown in Pig, 4,14. This transfer function could be reali- 
zed using digital filter technique and stored in memory 
of microcomputer. Taking the Fourier transform of digi- 
tized ISEMG and multiplying it with the transfer function, 
gives the Fourier transform of the force and its inverse 
Fourier transform gives muscular force. This was done 
with ISEMG picked up at different muscle contraction 
rates. ISEMG for muscle contraction rates of 4 N/sec3nd 
2 N/sec, are shown in Fig. 4.15 and Fig, 4.19, 2he 
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results of above stat-d procedure are shown in Fig. 4.15 to 
Fi9 • ^ ^ found that force determined using ISEMG 

(Fig. ^'^9* 4.22) is very similar to actual force 

recorded during exieriment (Pig, 4.23 and Fig, 4,24), 


4,5 ol S 2 3 3 < h 

The tran-.fer function of the black box found using 
both analo.; a:.i Ji.rital technique# can serve as a suitable 
device fo: mtd.suring muscular force from ISSMG picked up 
during locomotion. 
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CHA.PTER - 5 


SCOPE OF FURTHER WORK 


The basic objective of finding muscular force and 
ISEI^ correlation under non-isometric contraction of muscle 
is to determine muscular force by measuring SEMG. It has 
been explained earlier that during locomotion almost all 
muscle contractions are non-isometric. So, the transfer 
function model provided in the present work can serve as a 
suitable device for measuring force. 

Along with lifting and pushing and pulling, load 
carrying is a very common manual materials handling task. 

The metabolic energy requirsnents for different methods of 
load carrying have been studied under a variety of conditions. 
The general principle that has evolved from these studies is 
that in order to minimize energy expenditure the load should 
be located as closely as possible to the centre of mass of 
the subject. The total metabolic energy expenditure measure- 
ments provide no indication of local loading of particular 
muscles such as the muscles of the lower back. It is possible 
that a particular activity may require only a moderate incre- 
ase in total metabolic energy but cause fatigue in a local 
muscle group that has to produce high contractile level as part 
of accomplishing that activity. The amount of EMG output can 
be used as a measure of the relative tension in the muscle 
and from this, assessment could be made how different load 
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and position affect myogenic spinal joint compression and 
muscular fatigue. 

For the avoidarxre of low back problems investiga- 
tions can be done regarding proper lifting technique# body 
posture and maximum acceptable load of lifting. One of 
the possible contributing factors to injury is the effect 
of neuromuscular interaction in situations in which the 
individual has no knowledge of load magnitude prior to 
the actual lift attempt. Investigations can be done to 
determine the relationships among variables t^^pically associ- 
ated with lifting in order to minimize individual muscle 
forces and joint stresses. Studies can also be done to see 
whether past experience on specific task alters the muscular 
coordination and enhances overall movement skill by changing 
muscle action in a given situation. 
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